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COMPUTATIONAL HOLONOMY DECOMPOSITION OF 
TRANSFORMATION SEMIGROUPS 


ATTILA EGRI-NAGY1’2 AND CHRYSTOPHER L. NEHANIV^ 


Abstract. We present an understandable, efficient, and streamlined proof 
of the Holonomy Decomposition for finite transformation semigroups and au¬ 
tomata. This constructive proof closely follows the existing computational 
implementation. Its novelty lies in the strict separation of several different 
ideas appearing in the holonomy method. The steps of the proof and the 
constructions are illustrated with computed examples. 
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1. Introduction 

One of the fundamental concepts of science and computation is the notion of 
change: a system goes from a state to another state due to external manipulations 
or due to internal processes at various time-scales. If the set of states is a continuum 
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then we study continuous functions and thus we do analysis. If we have a set of 
discrete states then we do algebraic automata theory. A transformation semigroup 
{X, S) captures the concept of change in a rigorous and discrete way. It consists 
of a set of states X (analogous to phase space), and a set S of transformations of 
the state set, s : X X acting by x i-A x • s, that is closed under the associative 
operation of function composition. Writing S1S2 € S for the composite function 
Si G S' followed by S2 S S, we have x • (S1S2) = (x • Si) • S2, giving a (right) action of 
S on X. A fixed generating set for a transformation semigroup can be considered as 
a set of input symbols, therefore automata (without specifying initial and accepting 
states) and transformation semigroups are essentially the same concepts. 

Another fundamental technique of the scientific method is decomposition. The 
holonomy decomposition is a method for hnding the building blocks of a transfor¬ 
mation semigroup and compose them in a hierarchical structure. This composite 
semigroup has a structure that promotes understanding and it is capable of emu¬ 
lating the original transformation semigroup. Therefore, we say that the holonomy 
decomposition is a way of understanding transformation semigroups. 

Our aim here is to provide the simplest and most accessible proof for the holo¬ 
nomy decomposition theorem by giving a construction which is ‘isomorphic’ to its 
computational implementation [7, 5]. The novelty of this proof is the strict sep¬ 
aration of the several different ideas that appear in the holonomy decomposition. 
Both separating them from each other and from the technical details. 

1.1. General Ideas. There are four fundamental concepts used in the holonomy 
decomposition. First we state them in their generality to aid intuition, then give a 
short summary how they actually appear in the method. 

Approximation: gives less information about a system in a way that the 
partial description does not contradict the full description. 

Emulation: is a capability of one system producing the same dynamics as 
another one, not necessarily containing an exact copy. 

Compression: for repeated patterns stores the pattern once and record its 
occurrences. 

Hierarchy: is any system where the control information flows in one direction 
only and abstractions are natural operations. 

In the holonomy decomposition, we study the action on chains of increasingly 
smaller subsets of the state set, recovering the original transformations at the level 
of singleton subsets (approximation). Whenever the semigroup acts the same way 
on different subsets, we consider those subsets equivalent and only store the action 
on the equivalence class representatives (compression). These representative local 
actions are the building blocks and they are aligned according to a partial order 
(hierarchy). The chain semigroup and its encoded form, the cascade product can 
compute everything the original transformation semigroup can (emulation). 

1.2. Mathematical Preliminaries, Notation. A semigroup is a set S together 
with an associative binary operation S x S ^ S. A semigroup is a monoid if it 
contains the identity element. Let denote the monoid we get by adjoining an 
identity to S in case S is not a monoid. A transformation semigroup (A, S) is a 
finite nonempty set X (the state set) together with a set S of total transformations 
of X closed under function composition. The states are often denoted by a set of 
integers n = {1, ..., n}, and the transformations by the list of images [ji,j 2 , ■ ■ ■ ,jn]^ 
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where i i—>■ ji for i,ji £ n. The action a; i—>■ x • s on the points (states) x € X by 
transformations s € S naturally extends to set of points: P ■ s := {p ■ s | p S P}, 
P C X, s G S, and we have (P • si) ■ S 2 = P ■ (S 1 S 2 ), for si, S 2 G S. Similarly, the 
action can also be extended to sets of sets of points or to tuples or sequences of 
points or sets of points. 

The wreath product {X, S) I {Y, T) of transformation semigroups is the transfor¬ 
mation semigroup {X x Y,W) where 

W = {{sJ)\sGS,f GT^}, 
whose elements map X x Y to itself as follows 

{x,y) ■ (sj) = {x-s,y /(x)) 

toT X G X,y G Y. Here is the semigroup of all functions / from X to T (under 
pointwise multiplication). Note we have written y ■ /(x) for the element /(x) G T 
applied to y G Y. The wreath product construction is associative on the class of 
transformation semigroups (up to isomorphism) and can be iterated for any number 
of components. 

The size of the iterated wreath product grows rapidly by increasing the num¬ 
ber of components or by increasing their sizes. Explicit computation with wreath 
products is impractical. This motivates the definition of cascade products: effi¬ 
cient constructions of substructures of wreath products, induced by explicit de¬ 
pendency functions [6]. Essentially, cascade products are transformation semi¬ 
groups glued together by functions in a hierarchical tree. More precisely, let 
((Xi, Si), ..., (X„, Sn)) be a fixed list of transformation semigroups (here Si are 
semigroups and Xi the sets on which they act), and define dependency functions to 
be functions of the form 

di : Xi X ... X Xi-i — Si, for i G {2 ... n}. 

A transformation cascade is then defined to be an n-tuple of dependency functions 
d = (di,... ,dn), where di is a dependency function of level i. On the top level, di 
is simply an element of the semigroup Si. The transformation cascade d applied 
to (xi,..., x„) is defined coordinatewise by Xi ■ di{xi,..., Xi_i), applying the re¬ 
sults of the evaluated dependency functions, so that the cascade product can be 
regarded as a special transformation representation on the set Xi x ... x X„. The 
hierarchical structure allows us to conveniently distribute computation among the 
components, and perform abstractions and approximations of the system modelled 
as a cascade product. In the permutation group case it is basically the Schreier-Sims 
algorithm [14] put into product form [6]. 

1.3. Computational Tools. The constructive proof for the holonomy decompo¬ 
sition described here is implemented in the SgpDec [7, 5] software package for the 
Gap computer algebra system [10]. For the verification of the correctness of the 
software package we use a selection of transformation semigroups with interesting 
features and corner cases. We also have a shadow implementation of the algorithms 
based on partitioned binary relations in the kigen system [4]. 

1.4. Historical Notes. In Krohn-Rhodes theory, the holonomy method for cas¬ 
cade decomposition was originally developed by H. Paul Zeiger [21, 22], and sub¬ 
sequently improved by S. Eilenberg [8], and later by several others [3, 11, 13]. 
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Variants [2, 18], and generalizations of the theorem to the infinite case [9, 12] and 
to categories [ 20 ] were also studied. 

The term ‘holonomy’ is borrowed from differential geometry, since a roundtrip of 
composed bijective maps producing permutations is analogous to moving a vector 
via parallel transport along a smooth closed curve yielding change of the angle of 
the vector. 

The current proof is a prime example of the observation on the development of 
mathematics, that proofs turn into definitions (see the introduction of [19]), as the 
way we define the chain semigroup is the key argument of the previous proofs. 

2. Approximation 

For a transformation semigroup {X, S) we describe ways to approximate the 
states a; S V by subsets of X, and to approximate the transformations in S, the 
‘behaviour’ of the semigroup. 

2.1. Approximating states. What is the current state of the system? We can 
answer this question precisely by giving a single element, or we can give partial 
information by specifying a set of states with the condition that the current state 
is contained in the set. This way, any subset PCX such that x G P can be 
considered as an approximation of the state x. 

For a particular transformation semigroup we do not need to consider all such 
elements of the power set ViX), we can restrict to those that are generated by the 
semigroup action. 

Definition 2.1. The set Xs{X) = {X • s j s S S'} is the image set of the transfor¬ 
mation semigroup (A, S). 

Note that in general X itself and the singleton state sets are not necessarily 
included, so we may need to add them to the image set. 

Definition 2.2. The extended image set of the state set under the action of the 
semigroup is I's{X) = Xs{X) U {A} U {{xj [ x S A}. 

When approximating, we may be interested in doing it step-by-step. Since ap¬ 
proximations are subsets, we can build successive approximations by nested subset 
chains. 

Definition 2.3. A chain C is a subset of 'P(A) such that P G C and Q G C 
implies P C Q ot Q C P. A chain C is maximal if it is not properly contained in 
any other chain. We say that two chains C and D agree down to P P G C CT> 
and for all subsets Q with P C Q C X we have Q G C Q G T). 

Observation. Notice that S acts on subset chains in A: Since P C Q implies 
P ■ s C Q ■ s, necessarily C • s is a chain if C is. Moreover, (C • si) • S 2 = C • S 1 S 2 . 
However, the length of chains can become shorter under this action. 

As mentioned before, for the holonomy decomposition we do not need the full 
power set. However, we need the extended image set if we want to describe all 
necessary stages of approximating a state by maximal chains. 

Definition 2.4. Let C = C(A, S) denote the set of all maximal chains in I's(A). 

There is a surjective function 77 : C — » A mapping each maximal chain C to the 
element of its unique singleton {xj G C. We say C is a lift of x G A if r]{C) = x. 
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2.2. Approximating Transformations. A state x is lifted as a maximal subset 
chain starting from {a;}. Consequently, for lifting transformations we need to con¬ 
struct transformations mapping C to itself. However, simply acting on maximal 
chains, C i—> C • s is not a well-defined action on C, since C • s may not be maximal. 

Definition 2.5. A dominatinq chain of a chain C is a maximal chain D such that 

CCD. 

There can be more than one dominating chain. For instance, acting by a constant 
map on any chain would produce a singleton set, which can be dominated by all 
maximal chains containing that set. 

For any fixed s S S' we can define a (non-unique) mapping s : C —>■ C by 
s(C) = D, where D is any fixed maximal chain containing C • s. We can think 
of such an s as mapping the nested approximations C of t = 17 (C) to nested 
approximations s(C) of a; • s = ri(s{C)). We say s is consistent with chain structure 
if C and C' agree down to P then s(C) and s(C') agree down to P • s. 

One way to ensure this condition is to totally order I's(A), and for example 
choose its least member that can be included when building a dominating chain. 
We observe there is always at least one way to choose s so that it is consistent with 
chain structure. 

Lemma 2.6. If §1 and §2 mapping C to itself are consistent with chain structure, 
then so is the composite mapping S 1 .S 2 . 

Proof, ff maximal chains C and C' agree down to P G C n C' then, since Si is 
consistent with chain structure, si(C) and si(C') agree down to P • si. Since §2 is 
consistent too, we have that S 2 (si(C)) and S 2 (si(C')) agree down to (P • si) • S 2 = 
P-(siS 2 ). □ 

Definition 2.7 (Chain semigroup). Given a generating set A, for S, for each a G A 
we choose a consistent a and take S' = (a | a G A). Then we call the transformation 
semigroup (C, S) a chain semigroup. 

We say Oi • • • Ofc is a lift of s G S if s = Oi • • • for generators ai,... G A. 

By Lemma 2.6 it follows that any s = ai ■ ■ ■ aj- is consistent, i.e., 

Proposition 2 . 8 . All mappings in a chain semigroup (C, S) are consistent with 
chain structure. 

Remark 2.9. (1) We generally take just one lift a for each generator a of S to 

generate a chain semigroup, since one would often like S to be as small as 
possible. Different choices of lifts for the generators can result in different 
sized S. 

(2) Generally, there can be many different lifts s G S for fixed s in S, since 
s = ai ... Ok = a'l ••• a'^ does not imply di .. .au = 0 !^. . .a! g, although both 
are lifts of s. 

(3) There is a unique maximal chain semigroup obtained by taking all possible 
consistent s for s G S, and letting S be the semigroup they generate. 

In a sense chain semigroup contains approximations of (A, S). The rest of the 
holonomy decomposition is about putting an efficient notation (by embedding it 
into a wreath product) on this expanded semigroup. 
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3. Emulation 

We need to show that a chain semigroup emulates the original semigroup. 
Lemma 3.1. There is a surjective morphism of transformation semigroups 

{X,S)^{C,S). 

Proof. There is a semigroup homomorphism from S' to S' determined by a i-A a, 
where we recall that a is a generator of S. It is not hard to see this is well-defined. 
(And it follows, e.g., from Proposition 1.10 in [3]). Since 77(a(C)) = a; • a for 
X = r]{C), the action is respected. □ 

In the final form of the holonomy decomposition we will use the following notion 
of emulation. 

3.1. Division. One transformation semigroup divides another, (A, S) | (Y,T), if 
(A, S) is a homomorphic image of a substructure of (Y, T): precisely, there exists a 
subset Z (-Y and a subsemigroup U <T, with z ■ u € Z ioi dll z € Z,u G U, and 
a surjective function 9i : Z ^ Y and surjective homomorphism O 2 '■ U ^ S such 
that 6*1 (z • u) = 0i{z) ■ 92 (u) for all z & Z and u G U. 

4. Compression 

4.1. Equivalence of Subsets. On I's{X) we define an equivalence relation by 

P =s Q iff 3s, t G such that P = Q ■ s and Q = P -t. 

This is the equivalence relation of ’mutual reachability’ under the action of S, and 
the equivalence classes are the strongly connected components of (A, S) acting on 
I'siX). 

It is immediate that P =5 Q ==^ |P| = \Q\. As we will see S acts the 

same way on equivalent elements (see permutator and holonomy groups defined 
below), thus the equivalence classes provide the way to compress information in 
the decomposition. For each equivalence class there will be only one component in 
the hierarchical decomposition. 

4.2. Group Actions. For a subset P C X we have the stabilizer semigroup Sp = 
{s G S \ P ■ s = P}. If we restrict the action of the stabilizer to P we get the 
permutator group Gp. These groups are also called generalized Schiitzenberger 
groups [17]. 

In the holonomy decomposition we need the most coarse-grained approximation 
possible so we have to take another homomorphic image of Gp. Considering the 
inclusion relation (I's(A),C), we call a (lower) cover Pi of a non-singleton subset 
P G 1(A) a tile denoted by Pi -< P. The set of all tiles of P is denoted by T{P). 
These are the maximal subsets of P in I's(A). Obvious properties of tiles are: 

k 

P =[jPi, P,C Pj p^ = Pj 

i=l 

where Pi G T{P) and k = |T(P)|. Important to note that tiles of a set may overlap, 
so one should think of roof tiles as the analogy. 

The holonomy group Pip is the permutation group {P{P),Gp) made faithful. 
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Figure 1. An illustrative example of how moving around in an 
equivalence class induces permutations on the elements of the 
equivalence class. 


4.3. Constructing Holonomy Groups. If P =s Q, then there exist mappings 
mp^Q, niQ^p G S mapping P to Q bijectively {Q to P respectively), such that 
mp^QmQ^p is the identity map restricted to P and mq^pmp^Q is the identity 
restricted to Q (see e.g. [17]). 

It can be shown that if P =$ Q then Gp = Gq. Since there is a bijection between 
T(P) and T{Q), it follows that Hp = Hq. Moreover, ‘roundtrips’ of mappings in 
the equivalence class induce permutations on elements of the equivalence class (see 
schematic drawing on Figure 1). We can get the generators of Gp by contracting 
roundtrips of the form 


R 



P -^^ Q 


where P and Q are elements of the equivalence class of R and s is a generator of S 
mapping P bijectively to Q. 

5. Hierarchical Structure 

The output of the holonomy decomposition algorithm is a cascade product of 
transformation semigroups. So far we have established that the components of this 
cascade will arise from the holonomy groups of equivalence class representatives, 
but we still do not know how the components are put together in the cascade. 

5.1. Subduction. The inclusion relation is naturally defined on I'{X). In the 
subduction relation we also allow the sets to be moved by S. 

PQsQ 3s gS^ such that PCQ-s P,Q G T'{X), 

i.e. either P <G Q or we can transform Q to include P under the action of S. 
Therefore, subduction is a generalized inclusion, i.e. inclusion is subduction under 
the action of the trivial monoid. 
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It is easy to see that C 5 is a preorder: it is reflexive, since P Q P ■ 1, and it is 
transitive, since if P C Q • Si and Q C R ■ S2 then P C R ■ S2S1, thus P Cg R. 



Using a common technique for preorders, we define the =s equivalence rela¬ 
tion on I's{X) by taking subduction in both directions: P =s Q P C 5 

Q and Q Cg P. 

5.2. Positioning the components: Height and Depth of Sets. The height of 
a set Q e X'g{X) is given by the function h : X'g{X) —> N, which is defined by 
hs{Q) = 0 if Q is a singleton, and for |(5| > 1 , hg{Q) is defined by the length of 
the longest strict subduction chain(s) in the skeleton starting from a non-singleton 
set and ending in Q: 

hg{Q) = max((5i Cs ■ ■ ■ Cs Qi = Q), 

i 

where \Qi\ > 1. The height of {X,S) is h = hg{X). 

It is also useful to speak of depth values, which are derived from the height values: 

d{Q) = hg{X) - hg{Q) + 1. 

The top level is depth 1. 

Calculating the height values establishes the hierarchical levels in the decompo¬ 
sition, i.e. the number of coordinate positions in the holonomy decomposition is 
hs{X). 

Fact 5.1 (Depth never decreases). Let P € X'. Then d{P ■ s) > d{P). If d{P ■ s) = 
d{P) then P =g P ■ s. 

5.3. Positioned Chain Semigroup. By using the depth function, we can know 
align the members of those maximal chains on which the chain semigroup acts. 

Definition 5.2 (Positioned chain). For a maximal chain C 

X = PiD P 2 D P 3 D ...D Pk = {4. 

we take the associated positioned chain Cp°®. This is a vector of length hg{X) 
where the slots are empty (denoted by *) except that P^+i is in position d{Pi) for 
1 <i < k. For a positioned chain Cp°® the content at level i is 

This puts the members of chains into coordinate slots. By the maximality of the 
chain we have Cp°®[i] ^ Pi. Note that a positioned chain omits X, since it is not a 
tile of anything. 

We can identify the action of the chain semigroup with an action on positioned 
chains denoted by CP°®(X, S): 

Fact 5.3. {C,S) ^ {Cp°^,S). 

Proof. The positioned chains are in one-to-one correspondence with the maximal 
chains of C by the maps C Cp°®, since the only missing element of the chain in 
the positioned chain is X itself, so it can be added without any ambiguity when 
recovering the maximal chain. □ 
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At each level of depth we need to know how far the approximation proceeded so 
far, i.e. we need to know what subset of the state set are we acting on at the given 
depth. The value at the position is a tile, and tiles can belong to more than one 
set, so we need to look back to the first concrete value above. 


Definition 5.4 (state of approximation). 


a.(CP°^) 


a: if f = 1 

Cposjjj otherwise, where j 


maxj {CP°®[j] 7 ^ * and j < i} , 


1 < i < hs{X). 

Since ai only depends on where j < i, Ui is well-defined on prefixes of C of 

length at least i — 1. Moreover, we define a(CP°®) = (ai(CP°®),..., (Cp°®)). 

Lemma 5.5. For all maximal chains C and 1 <i < hs{X), 

CP°'^[i] G T{a,{CP°^)) 


when CP°®[z] ^ *. 

Proof. This is immediate from the definition of positioned chains (Def. 5.2). □ 


Lemma 5.6. For all maximal chains C it always holds the d(ai(CP°®)) > i. 

Proof. If f = 1 then ai(CP°®) = X and d{X) = 1 so the statement holds. If z > 1 
then assume the statement holds by induction hypothesis for i so we show that it 
follows for z -I- 1: 

Case 1: (i(ai(CP°®)) = z then by the definition of positioned chains Cp°®[z] is 
a tile of ai(CP°®), so at level z -I- 1 the value of Oi+i will be this tile, which 
is of depth at least z -I-1. 

Case 2: (i(a,(CP°")) > z then a,+i(CP°") = a,(CP°") but still fi(ai+i(CP°")) > 
i -I- 1. 

□ 


When lifting a transformation s, we only need to act when we are on the right 
level, i.e. d(ai(s(CP°®))) = z. The next lemma shows that the action of a lifted 
transformation respects approximation. 

Lemma 5.7. For a transformation s € S and a maximal chain C, we have for all 
coordinate levels i 

(a,(CP-)).sCa,((s(C))P-). 

Proof. Let Pi = ai(CP°^) and Qi = ch(('S(C))p°®). Pi = Qi = X, so the statement 
is true for z = 1. 

By induction hypothesis, the statement holds for levels down to and including z. 
Trivially, C P, and Qi+i C Qi. We show that P^+i • s C 

If Qi+l = Qi then the statement holds since P^+i • s C Pj • s C Qi. Otherwise, 
Qi+i C Qi and by the maximality of the chain Qi+i e TiQi)- 

Case 1: If P^ • s = Qi, then d{Pi) < d{Qi) as Pi cannot be deeper than Qi 
and d{Qi) = i since we are on the right level and d{Pi) > i always holds. 
Thus we have Pi =s Qi and d{Pi) = d{Qi) — i. Also, d{Pi^i) > z -f I, 
therefore Pi ^ Pi+i. Finally, Pi+i ■ s G s(Cp°®) ^ P^+i • s C Qi+i. 

Case 2: If Pi • s C Qi since s(C) is a maximal chain containing Pi • s and 
Qi+i is a tile of Qi, so Pi - s Q Qi+i, whence Pi+i • s C Qi+i- 

□ 
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5.4. Holonomy Cascade Semigroup. We build a cascade product of the holo- 
nomy groups of {X, S). First the components. Let i?i,..., be the representative 
sets of depth i. Then the ith component of the cascade product is defined as the 
transformation semigroup 

n, = i%,Tr,) = (r(i?i) u • • • u TiRk) u ,i<i< hs{x). 

The set of states are the set of tiles of the representative sets of depth i. These tile 
sets may overlap, thus we need to take the disjoint union. This causes no confusion 
since for each positioned chain we know the current state of approximation, hence 
we know which set of tiles we need to choose from. 

The transformations come from the holonomy groups of the representatives of 
depth i. How does Hi act on 77? If P lies in the jth set T{Rj) of the disjoint 
union then {hi ,..., € Hi acts on P by applying hj and it acts on * trivially. 

Recall that Hi augments the group Hi with all constant maps on 77- Since (77, Hi) 
is a well defined transformation semigroup for 1 < i < hs{X), we can form their 
wreath product. 

Definition 5.8. We call Hi I ■ ■ ■ I Hd — H{X,S) the holonomy wreath product 
semigroup of {X, S). 

In practice, we only want a substructure of this potentially huge wreath product, 
so we need to construct a cascade product by giving explicit dependency functions 
in the transformation cascades induced by the generators of S. The maps are 
s i-A s I—)■ enc(s), where the final encoding describes s in terms of the corresponding 
representative set. 

5.4.1. Encoding and decoding. We encode the elements of a positioned chain, that 
are tiles of the current state of approximation, as tiles of the representative set of 
the corresponding height. If Cp°®[7] ^ * then 

enc(CP°'^)[j] = CP°"[j] • rup^p where P = aj(CP°'*), 

otherwise the encoded value is *. Since a is not recursive, encoding can also be 
done independently for any level. 

Decoding does the opposite, however we need to calculate the current unencoded 
state of approximation, therefore it is a recursive calculation. Let V = enc(CP°®), 
the tuple of coordinate values. If Vp°®[7] 7 ^ * then 

dec(VP°®)[j] = VP°®[j] • rrip^p where P = aj(dec(VP°®)), 

otherwise the encoded value is *. These are bijective maps, thus dec(enc(CP°®)) = 
CP°" and enc(dec(V)) = V. 

5.4.2. Dependency functions. For s in a chain semigroup S, let’s define enc(s) to 
be the transformation cascade given by the dependency functions 

(enc(s)) ■ : Ti X • • • X 77-i ^ H,. 

Let’s fix a positioned chain Cp°®, and thus P = ai(CP°®), Q = ai(s(CP°®)) and 
V = enc(CP°’*). By Lemma 5.7, these state approximations satisfy P ■ s C Q. 

We need to define the value of the dependency function (enc(s)) ■ on (Hi,..., H-i), 
the prefix of H: It is constant * G Hi unless we are on the right level, i.e. i = d{Q), 
in which case we have a constant map (reset) to a tile or a permutation. 

Precisely, if z ^ d{Q), then let enc(s)i(Vi,..., H-i) = constant * G Hi. 

There there are two possibilities when i = d{Q): 
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Permutation: If the chain action satisfies P ■ s = Q. The encoding of s at 
depth i on chains that agree with Cp°® up to depth i — 1 is rrip^p s tuq^q, 
a permutation oi P = Q, therefore 

enc(s)i(Vi,...,yi_i) = mp^pSiriQ^Q, 

and this is in the component of the holonomy group of P by the definition 
of holonomy groups, and has identities elsewhere according to the disjoint 
union action. 

Reset: If the chain action satisfies P ■ s C Q according to s we take the tile 
enc(s(CP°®))[i] of the representative Q and let 

enc(s)i(pL,... = constant(enc(s(CP°®))[i]). 

Since s is consistent with chain structure this constant is the same for 
all chains that agree with Cp°® up to depth i — 1. Again, the value of the 
dependency function is in Hi by the definition of the holonomy permutation- 
reset transformation semigroups. 

It is clear that enc(s)i is well-defined since any s in the chain semigroup is 
consistent with chain structure, and enc and dec are defined level-by-level on chains 
(same prefix gives same result). Therefore we have an element enc(s) of the wreath 
product, i.e., enc(s) £ 'H(X, S) 

Theorem 5.9. (C, S) = (CP°®, S) ^ Hii ■ ■ ■ i Hd = P-iX, S), where d = hs{X). 

The image of such an embedding is called a holonomy (decompostion) cascade prod¬ 
uct. 

Proof. The isomorphism was shown in Fact 5.3. We show enc is an embedding of 
transformation semigroups from (Cp°®, 5') to the wreath product. For the states, 
enc(CP°®) C 7i X • • • X 7d holds trivially. We need to show that if V = enc(CP°®) £ 
enc(CP°®) then enc(s)(V) = enc(s(dec(V))). 

By looking at the ith position for each 1 < i < d, it i ^ d{s{ai{CP°^))) then 
s(CP°'^) cannot have a tile in position i, it follows that * = enc(s(CP°'’))j, which is 
equal to 

Vi ■ enc(s)i(Vi,..., Vi_i) = Vi ■ constant *, 

as required. 

Otherwise, i = ci(s(ai(CP°®))), and we have two cases. If enc(s)i(Fi, ■ • • , Vi-\) is 
a constant map to a tile, then the definition of enc(.§) yields 

Vi ■ enc(s)i(Fi,... ,Vi-i) = Vi ■ constant(enc(s(CP°®))[z]) = enc(s(CP°®))[i], 
as required. Otherwise, the component action is a permutation, and then 
enc(s)(V)[*] = Vi ■ enc(s)*(Fi,..., Vi-i) 

= enc(CP°")[t] • mp^pSmQ^Q 
= CP°’*[i] • mp^p m-p^p s tuq^q 
= CP°®[f] • suiq^q 
= enc(s(CP°®))[f] = enc(s(dec(V)))[z] 

by the property that rup^-p m-p^p = Ip, the identity map on P, hence on its set 
of tiles, where P = ai(CP°®) and Q = Q;i(s(CP°®)). 
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Since this holds for all 1 < i < d, we have 

enc(s)(enc(CP°®)) = enc(s(CP°®)). 

It follows that 

(enc(s 2 ) o enc(s"i))(enc(CP°®)) = enc(s 2 )(enc(si)(enc(CP°’*)) 

= enc(s2)(enc(si(CP°®))) 

= enc(s2(si(CP°^))) 

= enc(s2 o si)(enc(CP°®)). 

Thus, enc is clearly an (injective) semigroup homomorphism. Whence, (enc(CP°®), enc(5')) 
is a isomorphic to (Cp°®, S). □ 

By Lemma 3.1 and Theorem 5.9, we have 

Corollary 5.10 (Holonomy Decomposition Theorem). A finite transformation 
semigroup {X, S) divides its holonomy wreath product 

{x,s) \nii---ind = n{x,s), 

where d = hs{X). 


6. Computational Complexity 

The holonomy decomposition algorithm given here enumerates the image set 
Is{X) of the state set X. The worst case is enumerating the powerset with 
elements. It is easy to conclude that the algorithm given has time complexity 
at least exponential in the number of states (cf. Maler [18]). Moreover, by the 
Krohn-Rhodes prime decomposition theorem [15, 16], every simple group divisor of 
a hnite semigroup must occur as a divisor of any cascade decomposition. Therefore 
it follows that a finite automata has no nontrivial subgroups (i.e., is aperiodic) if and 
only if all its holonomy groups are trivial. The results of Cho and Huynh [1] show 
that aperiodicity is PSPACE-complete, so it follows immediately that computing 
the holonomy decomposition is PSPACE-hard. 

In practice we can calculate with huge semigroups (of size hundreds of thousands 
of elements). The size of the state set and the size of generator set or of the 
semigroup do not necessarily give a good guide to computational complexity in 
practice. It would be interesting to find the appropriate features and parameters 
and do parametrized complexity analysis for holonomy decompositions. 

7. Computed Examples 

Example 7.1. As a minimalistic but non-trivial example, let (3,5) be the trans¬ 
formation semigroup generated by si = [2,1,3] and S 2 = [1,2,2]. From Figure 2 
we can read off the maximal chains: {{1, 2,3}, {1, 2}, {1}}, {{1, 2,3}, {1, 2}, {2}}, 
{{1,2,3}, {3}}. Let’s see how from t = S2S1 = [211] we construct t acting on 
the chain representing state 1, i.e. doing the action on the members of the chain, 
removing duplicates then finding a dominating chain. 


chain C 

C-t 

C-t 

{1,2,3} 

{1,2} 

{1,2,3} 

{1,2} 

{1,2} 

{1,2} 

{1} 

{2} 

{2} 
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Figure 2. The skeleton of a semigroup acting on 3 points (Ex¬ 
ample 7.1). The nodes are the elements of I's{X). The boxes 
are the equivalence classes, the rectangular nodes the chosen rep¬ 
resentatives of a class. Shaded equivalence classes have nontrivial 
holonomy groups. The arrows point to the tiles of a representative 
set, the labels denote sequences of generators taking the set to its 
tile. Dotted arrows indicate tiles that are not images. On the side 
depth values are indicated. 


In this very small example we have only a single dominating chain. 


Example 7.2. Let (6,5') be the transformation semigroup generated by transfor¬ 
mations Si, ..., Sg- 


creates the image {1, 2, 3}, 

is the transposition (4,5) and gives the image {4,5,6}, 
is a cycle on (4, 5,6}, 
creates the image (4, 5}, 
maps {4,5,6} to {1,2,3}, and {1,2,3} to {4}, 
is a cycle on {1, 2,3}, 
and its basic properties are |5| = 138, and I's{X) = 19. The ‘skeleton’ of its 
holonomy decomposition is depicted on Figure 3. 

For C = {{1,2,3,4, 5,6}, {1,2,3,4}, {2,4}, {2}} we have 


51 = [1 2 3 1 1 1] 

52 = [4 4 4 5 4 6] 

53 = [4 4 4 5 6 4] 

54 = [4 4 4 4 5 5] 

55 = [4 4 4 1 2 3] 

56 = [2 3 1 4 4 4] 


depth 

^pos 

enc(CP°") 

a(CP°") 

1 

{1,2,3,4} 

{1,2,3,4} 

{1,2,3,4,5,6} 

2 

{2,4} 

{2,4} 

{1,2,3,4} 

3 

{2} 

{1} 

{2,4} 

4 

* 

* 

{2} 

5 

* 

* 

{2} 


demonstrating that an encoded positioned chain is not necessarily a chain. 


Example 7.3. The full transformation semigroup Ta has a canonical generating 
set consisting of two permutations (transposition and cycle) and an elementary 
collapsing. Figure 4 shows how these generators act on the set of maximal chains. 
The generators [2,1,3] and [2,3,1] are permutations of {1, 2, 3} that map maximal 
chains to maximal chains. The lifts of these transformations to the chain semigroup 
are thus exactly as shown the Figure 4 and hence unique. The transformation 
[1,2,1] gives subsets chains {{1, 2}, {1}} and {{1, 2}, {2}} that miss the full state set 
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Figure 3. The skeleton of a semigroup acting on 6 points (Exam¬ 
ple 7.2). Interesting feature is that (i({l,4}) < d({l,2,3}). 


itself. Thus, [1, 2,1] does not map maximal chains to maximal chains, but it maps 
{1,2,3} to {1,2}. For any lift of [1,2,1]: All maximal chains trivially agree down 
to {1, 2, 3} so it must map to a chains agreeing down to {1,2} = {1, 2,3} • [1, 2,1], 
thus the lift the of [1,2,1] maps each maximal chain C to {1,2,3} D {1,2} D 
{ 77 (C) • [1, 2,1]}, and so is uniquely determined. 

However, having a unique dominating chain or unique lift is not a general prop¬ 
erty. Constant map c = [3,3,3] produces the chain {{3}} for which any maximal 
chain containing {3} is a dominating chain. Since any two maximal chains Ci and 
C 2 both start with the top set X = {1, 2, 3}, they agree at X and so, by consistency 
c(Ci) and c(C 2 ) must agree down to A • c = {3}. That is, c(Ci) = c(C 2 ), and 
c is itself a constant map. Here there are two choices, {1,2,3} D {1,3} D {3} or 
{1, 2, 3} D {2, 3} D {3}, for the constant value of c. 

The same argument applies to lifting any constant map in this holonomy method: 
the lift of a constant to the chain semigroup yields a (non-unique) constant. 
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c 

C-[2,l,3] 

{1,2,3} 

{1,2,3} 

X 1 \ 

X 1 \ 

{1,2} {1,3} {2,3} 

{1,2} {2,3} {1,3} 

/\ /\ /\ 

/\ /\ /\ 

{1}{2} {1}{3} {2} {3} 

{2}{1} {2} {3} {1}{3} 

C-[2,3,l] 

C-[1,2,1] 

{1,2,3} 

{1,2} 

/ 1 \ 

X 1 \ 

{2,3} {1,2} {1,3} 

{1,2} {1} {1,2} 

/\ /\ /\ 

/\ /\ /\ 

{2} {3} {2}{1} {3}{1} 

{1}{2} {1}{1} {2}{1} 


Figure 4. Action of the canonical generators (transposition, cy¬ 
cle, elementary collapsing) on the maximal chains. 
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Appendix A. Notation 


Cs 

C, D 

CPOS, CP°"[i] 

Cp, Cx, c 

{cA) 

QpOS 

Gp 

Ez 

Hp 

hs{P),d{P) 

’^P-s-P 

P,QeIsJ^X),l'siX) 

P=sQ,P 

Sp 

r{Q) 

V 

(A, 5), (y,r) 

x,y,z G X 


subduction relation 
tile of relation 

current state of approximation at depth i for a positioned chain 
chains 

positioned chains, content of position i 

maximal chains from P, {a;} to X, all maximal chains. 

chain semigroup 

all positioned tile chains 

the permutator (generalized Schiitzenberger) group of P 

holonomy permutation group of P, {T{P),Gp) made faithful 

holonomy permutation-reset transformation semigroup of P 

height, depth of a set 

mapping from and to a representative 

images, image set, extended image set 

equivalence, representative element 

setwise stabilizer semigroup of P 

the tiles of Q 

encoded coordinate values (tiles of representatives) 
transformation semigroups 
states, stat set 
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